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The heat transfer performance of inhomogeneous parallel plate heat exchangers in transient operation is 
investigated using an established model. A performance parameter, denoted the Nusselt-scaling factor, is 
used as benchmark and calculated using a well-established single blow technique. A sample of 50 ran¬ 
dom stacks having equal average channel thicknesses with 20 channels each are used to provide a statis¬ 
tical base. The standard deviation of the stacks is varied as are the flow rate (Reynolds number) and the 
thermal conductivity of the solid heat exchanger material. 

It is found that the heat transfer performance of inhomogeneous stacks of parallel plates may be 
reduced significantly due to the maldistribution of the fluid flow compared to the ideal homogeneous 
case. The individual channels experience different flow velocities and this further induces an inter-chan¬ 
nel thermal cross talk. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

MicroChannel heat exchangers show promise for applications 
that require high cooling power density. It has been suggested that 
liquid cooling using a microchannel heat exchanger is the best sui¬ 
ted technique for cooling electronics as transistor density contin¬ 
ues to increase [1], Microchannels have also been applied to a 
wide variety of applications, such as cryocoolers, dehumidifiers, 
Stirling engines, solar power, electronics cooling and magnetic 
refrigeration [2-8]. 

Despite the high theoretical performance, microchannel heat 
exchanger performance in actual devices has often been reported 
in experiments to be lower than expected [9]. Several explanations 
for the relatively large deviations observed have been suggested, 
such as whether the fluid continuum assumption breaks down, 
the influence of surface roughness in the channels etc. In Ref. [9] 
these issues are reviewed and it is concluded that for incompress¬ 
ible, single phase laminar flows with aqueous fluids no new phys¬ 
ical phenomena occur in microchannels. This is supported by 
careful experiments performed on single-channel tubes and square 
channel heat exchangers in the microchannel range [10,11]. 

It was shown that flow maldistribution in fluid manifolds can 
reduce the microchannel heat transfer performance in a micro- 
channel heat exchanger [12], as channels near the edges of the heat 
exchanger do not receive as much fluid flow as those in the center. 
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The flow and temperature distributions in two parallel microchan¬ 
nels were studied numerically when obstructions such as bubbles 
or debris were placed in one of the channels [13]. The outlet tem¬ 
perature profile was shown to be affected by an obstruction in one 
of the channels, but the accompanying change in heat transfer per¬ 
formance was not quantified. 

Another explanation of the observed heat transfer degradation 
is non-uniform plate spacing [14] in the microchannel stack. As 
heat transfer performance is inversely proportional to plate spac¬ 
ing, the plate spacing is reduced as much as possible for many 
applications. As dimensions become smaller the relative manufac¬ 
turing tolerance can become significant. It was shown that varia¬ 
tions in fluid channel heights can dramatically reduce the 
effective heat transfer in a heat exchanger consisting of a stack of 
microchannels because larger fluid channels receive a dispropor¬ 
tionately high fluid flow while smaller channels are starved for 
flow [14J. The effect was demonstrated through numerical model¬ 
ing as well as experiments on passive regenerators. It was also 
shown that thermal cross talk, heat transfer between adjacent fluid 
channels through channel walls, equalizes temperature differences 
between fluid channels caused by differences in channels heights. 
Therefore, cross talk between channels generally reduces the neg¬ 
ative impact on heat exchanger performance due to a non-uniform 
distribution of fluid channels. 

This article studies the effects of thermal cross talk in more de¬ 
tail and shows its effects on microchannel heat exchanger perfor¬ 
mance under transient operation. It will be shown that the effect 
of the cross talk can vary greatly with geometry, material 
properties, and operating conditions. We show that a non-uniform 
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microchannel stack can fall into three categories from a thermal 
modeling standpoint: isolated channels, cross talk dominated, 
and effectively uniform. In the case of isolated channels, the ther¬ 
mal interaction between the channels is so low in relation to other 
phenomena that the channels can be modeled as individual iso¬ 
lated channels operating in parallel with an acceptable accuracy. 
In cross talk dominated operation, a detailed model of a stack of 
non-uniform flow channels is needed in order to obtain accurate 
results. In the last category the cross talk is sufficiently high in rela¬ 
tion to other heat transport, the heat exchanger can be assumed to 
be a uniform stack. It is shown that a given microchannel stack can 
fall into any of these three categories of cross talk depending on the 
fluid flow rate or thermal conductivity of the solid. A detailed 
numerical model of a microchannel stack is used to calculate the 
temperature response when the heat exchanger is subjected to a 
step change in temperature at its inlet. A technique developed 
for characterizing heat transfer in porous beds [15] is used to cal¬ 
culate the bulk heat transfer in the entire regenerator and compare 
it to an ideal microchannel stack with uniform fluid channel 
heights. All results generated here use a transient response tech¬ 
nique that was developed primarily for regenerators rather than 
steady state heat exchangers such as heat sinks. However, the gen¬ 
eral conclusions should also apply to heat exchangers in general. 

2. Heat transfer in inhomogeneous parallel plate stacks 

The general problem of determining the heat transfer character¬ 
istics of a heat exchanger, here specifically parallel plate stacks, in¬ 
volves thermal conduction in a solid material (of which the 
parallel plates are made), thermal conduction and convection in a 
heat transfer fluid in intimate contact with the solid and the heat 
transfer between the two media, in the idealized case the stack is as¬ 
sumed perfectly homogeneous, and it is therefore sufficient to con¬ 
sider half a solid plate and half a fluid channel due to symmetry. 
Numerical models of such systems have been reported including 
one-dimensional models where the direction of the flow (or axial 
direction) is the spatial dimension resolved and, in certain cases, 
heat transfer transverse to this direction is accounted for via a 
Biot-Fourier number approach [ 16], In other, 2-dimensional, models 
both the axial direction and the direction transverse to the flow (the 
y-direction in Fig. 1) are numerically resolved for the purpose of 
studying, e.g., low-thermal conductivity materials [17,18]. 

It is clear that if the stack is inhomogeneous, e.g., due to varying 
channel thicknesses and/or plate thicknesses that the symmetries 
just described will not apply. Each fluid channel will have a differ¬ 
ent mean fluid velocity, which will cause varying heat transfer in 
the transverse direction between adjacent channels and plates. 
This effect is referred to as thermal cross talk. It is non-trivial to 


provide a general estimate of the influence of this cross talk and 
in the following Section 2.1.3 a term that quantifies the signifi¬ 
cance of cross talk for a given operating condition is proposed. 

The other influence that inhomogeneity in a parallel plate stack 
has on the heat transfer characteristics is, of course, due to the vary¬ 
ing heat transfer conditions in each individual channel caused by the 
varying flow velocities (or flow maldistribution) and thus varying 
convection coefficients. Larger channels will generally have poorer 
heat transfer than smaller ones and at the same time significantly 
larger fractions of the total flow rates thus causing the heat exchan¬ 
ger to perform less efficiently than it would had it been a perfect 
stack with channel thickness equal to the average spacing of the 
inhomogeneous stack [ 14[. These two effects are investigated apply¬ 
ing the numerical model described in Section 2.1.2 below. 

2.3. Governing equations 

2.3.3. 3D homogeneous regenerator equations 

Porous medium regenerators are often modeled as a solid phase 
interacting with a fluid phase with uniform flow. The geometry of 
the regenerator is characterized by correlations for pressure drop, 
heat transfer and other correlations rather than by solving detailed 
conduction and fluid flow equations to determine them directly. 
The one-dimensional partial differential equations for the solid 
and fluid are given below. For the fluid the equation is 

m { c f ^ + hA m (T f - T s ) + p f A c ec f ^ - k f A c 0 = 0 (1) 

where T is temperature, p is density, c is specific heat, k is thermal 
conductivity, h, is the heat transfer coefficient, A H t is the area for heat 
transfer per unit length, £ = A is the porosity, rh f is the fluid mass 
flow rate, t is time and A c is the cross sectional area. Subscripts s and f 
indicate solid and fluid, respectively. The x-direction is defined as the 
direction of the flow (the axial direction; see Fig. 1 ). The terms repre¬ 
sent (in order from left to right) the enthalpy change of the flow, heat 
transfer from the fluid to the solid, energy storage, and axial conduc¬ 
tion in the fluid. Viscous dissipation due to pumping losses is ignored. 
Eq. (1) cannot directly account for in-homogeneous regenerator 
geometries because the fluid flow would be a function of they direc¬ 
tion. However, the correlations used to determine heat transfer, pres¬ 
sure drop etc. can be modified to account for in-homogeneous effects. 
The 1D solid material governing equation is 

M H t(7> - T s ) + M c 0 = (1 - e)AcP s Cs ^ (2) 

The terms represent heat transfer from the fluid to the regenerator, 
static axial conduction and energy storage. 



Fig. 1 . Schematic of the parallel plate stack geometry: (a) 3D perspective of the stack, (b) the geometry modeled in the detailed model. Each plate has the same thickness 
whereas the flow channels may have varying thicknesses. The x-direction is the direction of the flow and the y-direction is referred to as the transverse direction. 
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The two equations for a homogeneous regenerator, (1) and (2) 
can be non-dimensionalized by introducing the following terms 



(3) 


The characteristic time is z and the length of the heat exchanger, L, 
is used as a characteristic length scale. The hot and cold tempera¬ 
tures are denoted T H and T c , respectively. 

Substituting the terms in Eq. (3) into Eqs. (1) and (2) and rear¬ 
ranging give the non-dimensionalized fluid equation 


dx* 


+ NTU(d f - 8 S ) + 


9? 90f 
(p dt* 


1 


d 2 6f 


Pe dx* 2 


= 0 


(4) 


Combining Eqs. (1) and (2) and non-dimensionalizing them gives 


90s 

at* 


^ +9i 


ar 


a 2 ^ q> ctdf 
dx* 2 + Pe dx* 2 ' 


(5) 


The number of transfer units in the regenerator is defined as 

NTU = Mh^ (6) 

mfCf 


and the Peclet number is defined as 


Pe = Pf^ 

fef 


(7) 


where u is the average velocity of the fluid. The utilization is defined 
as 


nifCft 

CP ~ p s A c Lc s (\ - £) ’ 
where the mass flow rate is 


( 8 ) 


rhf = uep f A c . (9) 

The ratio of thermal capacity of the fluid and solid domains is 


91 = 


P f Cf£ 


p s c s (l - e) 

Finally, the Fourier number is 
K t 


Fo = 


L P s c :S (1 - e) 


( 10 ) 


( 11 ) 


2.J.2. Inhomogeneous stack model 

A parallel plate heat exchanger is, in the ideal case, straight-for¬ 
ward to analyze accurately in terms of its heat transfer character¬ 
istics. Flowever, in the case where the fluid channel thicknesses 
vary the coupled heat transfer problem of the solid and the fluid 
may become quite involved. The driving force of this is the asym¬ 
metry created by the distribution of channel thicknesses, which 
generally results in varying fluid velocities in the individual chan¬ 
nels- assuming the pressure drop across the heat exchanger is the 
same for each channel. Larger channels will dominate smaller 
channels due to their higher flow velocity. This may be seen from 
the analytical expression for the flow velocities [14]: 


; = V 




( 12 ) 


where u, is the mean fluid velocity of the ith channel, V' is the total 
volumetric flow rate of the stack and N is the number of channels in 
the stack. 


The number of transfer units describing the heat transfer be¬ 
tween solid and fluid includes the convective heat transfer coeffi¬ 
cient, h. This is defined through 



(13) 


where Nu is the Nusselt-number and d h is the hydraulic diameter of 
the heat exchanger. For the present case, flat plates and incompress¬ 
ible heat transfer fluid, the Nusselt number is virtually constant and 
has a value close to 7.54 [22]. The hydraulic diameter is d h = 2H f 
with H f denoting the fluid channel thickness. Due to varying chan¬ 
nel thicknesses, h will vary from channel to channel. 

A larger fluid channel has a lower heat transfer coefficient (Eq. 
(13)) combined with an increased flow velocity that results in a 
degradation of the heat exchanger performance. The thermal cross 
talk may help to restore some of the heat transfer efficiency of the 
heat exchanger [14]. 

In order to capture these different effects, it is necessary to ap¬ 
ply a detailed numerical model. Such a model was described in 
Refs. [23,14] for parallel plate heat exchangers with varying flow 
channel thicknesses. The geometry modeled is shown in Fig. 1. 
Thermal conduction is solved for in both the x- and y-directions 
and the fluid convection is included in the flow channel domains. 
The coupled equations solved for are 


fdT f dT f \ , 

PfCf V at + Ux ox ) = kfS/ Tf 
dT s . v7 2'r 
PsCs "aF = fcsV Ts ' 


(14) 

(15) 


The flow velocity profile is assumed to be fully developed at the 
channel entrance with a parabolic profile. The model is run with 
boundary conditions matching a single-blow, which means that 
the solid and fluid are initially at the same temperature, T init and 
the inlet fluid has a temperature T in . In this way the temperature 
breakthrough curve at the outlet of the stack is recorded as a func¬ 
tion of time. An example of such a curve is given in Fig. 2(a). 

This curve, T ou t(t), may be used to analyze the quality of the heat 
exchanger in terms of its efficiency [15]. Several parameters may 
be derived from the breakthrough curve that can be used to char¬ 
acterize the heat transfer performance. Here, the time difference 
between 20% and 80% of the breakthrough curve are used (see 
Fig. 2(b) and Refs. [14,15] for more detail). 

As described in detail in Ref. [14] a single-channel model is used 
to evaluate the heat transfer performance of the inhomogeneous 
stack. In the single-channel model it is possible to vary the thermal 
contact between the solid and fluid through a variable contact 
resistance, R. In this way the heat transfer between solid and fluid 
may be reduced and that leads to a Nusselt-number scaling factor 
[14]: 


Nu sc , = 


1/h 

1 /h + R' 


(16) 


Fig. 2(b) gives the characteristic breakthrough time (S) and the 
steepest gradient (M) as a function of the Nusselt-scaling factor 
found at a given flow rate from the single-channel model. In this pa¬ 
per, the detailed model is used to determine the time difference be¬ 
tween the 20% and 80% breakthrough for stacks of varying solid 
conductivity and channel thickness variation and match those val¬ 
ues to an equivalent single channel, or uniformly stacked, heat ex¬ 
changer model to determine the effective bulk heat transfer 
coefficient for the entire stack. 

It follows from Eqs. (6) and (16) that the NTU for an inhomoge¬ 
neous stack at a given geometry and flow rate is influenced by the 
derived Nusselt-scaling factor resulting in an NTU for an inhomo¬ 
geneous stack: 
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(a) (b) 


Fig. 2. Example of the analysis scheme applied for evaluating the heat transfer performance, (a) Shows the temperature breakthrough curve as a function of time. The 
steepets gradient and the time between 20% and 80% are indicated in the figure, (b) The S (breakthrough time) and M (steepest gradient) values as a function of Nusselt scaling 
factor (defined in Eq. (16)) for a fixed flow rate derived from the single-channel model. 


Table 1 

Central numerical properties/settings used in the simulations. 


Description/property 

Value 

Number of gridpoints in the x-direction (n x ) 

40 

Number of gridpoints in the y-direction in each 
domain (n y ) 

10 

Maximum timestep (At max ) 

0.01 s 

Numerical software 

Comsol Multiphysics 

4.2 

Solver package 

PARD1SO 


can be investigated by examining the heat transfer between chan¬ 
nels with unequal fluid velocities in channel i as 

Qd = Tf ' + n ~ Tf ' ( 18 ) 

Kcr 

where Kct is the thermal resistance between two fluid channels. The 
thermal resistance between the two fluid streams will include a 
convection resistance term, but for this work, R CT is assumed to 
be the thermal resistance associated with steady state conduction 
through the plate: 


NTU stack = 


Nu sc \NukfA m L 

d h rh f Cf 


(17) 


H s 

2A m Lk s 


(19) 


By reducing the inhomogeneous stack down to a homogeneous 
stack with a reduced NTU value, the complex geometry of non- 
homogeneous heat exchanger can be modeled in a ID numerical 
model, greatly increasing computational efficiency while preserving 
accuracy. Details on the numerical implementation may be found in 
Ref. [14]. Table 1, however, provides an overview of the important 
numerical parameters. 

2A3. Thermal cross talk between fluid channels 

The non-dimensionalization presented in Eqs. (3)—(11) is only 
valid for the idealized case where the heat exchanger is homoge¬ 
neous and temperature variations transverse to the flow direction 
can be ignored. 

In a homogeneous plate stack, the velocity in each channel is 
equal and heat transfer from the center of each plate to the fluid 
is symmetric about its centerline as shown in Eq. (12). When un¬ 
equal fluid velocities exist at either side of a plate, the heat transfer 
process becomes more complex. Rather than having equal convec¬ 
tion at both sides of the plate, there will be higher convection at 
one side and the high velocity fluid channel can transfer heat 
across the plate to the low velocity channel under transient oper¬ 
ation. This cross talk process is a complex thermal process that de¬ 
pends on transport properties of the solid, convection at the plate 
surface, and the time. This means that the fully detailed cross talk 
model would depend on Fo and other parameters as suggested by 
Eq. (5), and would probably have to be evaluated for a specific 
operating condition. Characterizing the cross talk in a complete 
manner would require an analysis such as in Ref. [16] and is out¬ 
side the scope of the present work. However, the general behavior 


where H s denotes the plate thickness. If a cross talk term were 
introduced into Eq. (1) and then non-dimensionalized, the resulting 
term would be in a form of a number of transfer units cross talk: 


NTUct 


1 

RcrrhfCf' 


( 20 ) 


When Rcr decreases the thermal cross talk becomes more important 
and may help to reduce the negative impact on the performance 
due to flow maldistribution. This happens if the solid thermal con¬ 
ductivity is increased and/or the plate thickness is decreased (see 
Eq. (19)). At large values of Rcr the cross talk becomes negligible 
and each flow channel may be treated as an individual channel, 
which makes analysis simpler [19-21]. 


2.2. Parameter variation 


It is clear that not only does the fact that the fluid channels have 
different thicknesses cause a degradation of the heat transfer per¬ 
formance due to varying flow velocities, but the degree of variation 
between the individual channel thicknesses may also be very 
important and the distribution within a given stack likewise. It 
is, therefore, in order to state anything general (given a certain 
parameter space in terms of geometry and transport properties) 
necessary to employ a certain amount of statistics to the analysis. 

In real applications the distribution of the channels will gener¬ 
ally be random with some overall standard deviation from the 
mean of the stack. Here, a total of 50 stacks each having 20 chan¬ 
nels (and thus also plates) were chosen to provide the base for the 
analysis. Each stack has the same average channel thickness but 
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Table 2 

Fixed parameters used in the detailed model. 


Parameter 

Cf (J/kg K) 

(J/kg K) 

Hf (Pa s) 

Pf (kg/m 3 ) 

Ps (kg/m- 3 ) 

k f (W/m K) 

L (m) 

Hf mean (mm) 

H s (mm) 

Value 

4200 

300 

0.001 

1000 

7900 

0.6 

0.04 

0.2 

0.4 


Table 3 


Parameters varied in the detailed model. 

Parameter 

V'/N [10- 6 m 2 s'] 

<r(*l 

k s (W/(m*K)) 

Value range 

2.5-60 

5-25 

0.5-400 


the distribution is different for each of the 50 stacks. The individual 
distributions for each stack were drawn randomly from a normal 
distribution and the same 50 distributions applied in all cases. 
The 50 distributions were scaled to keep their mean value but alter 
their standard deviation, a. 

From each stack, given a certain set of operating parameters, the 
respective S values were derived and the mean value of these 
found. This final mean S value was then used to derive the Nusselt 
scaling factor (Eq. (16) and Fig. 2(b)). It was found that using be¬ 
tween 40 and 50 stacks gives consistent results whereas only using 
20-30 stacks introduces some amount of noise to the results due to 
an insufficient statistical base. 

Clearly, heat exchanger stacks comprise a vast parameter space 
and only a selected part of this is treated here. Parameters that are 
kept fixed are given in Table 2 and those that are varied are given 
in Table 3. 

The choice of parameters is based on a heat transfer fluid sim¬ 
ilar to water and a solid material similar to gadolinium. The prod¬ 
uct of the mass density and specific heat is very similar when 
considering metals like gadolinium, aluminum, brass and copper 
(in the range 2.4-3.4J/g I<). These are common metals used in var¬ 
ious heat exchanger applications. Their range of thermal conduc¬ 
tivity is somewhat broader and that has been reflected in the 
choice of the variation of fc s ; see Table 3. The variation of the volu¬ 
metric flow rate, V', is chosen to give Reynolds numbers ranging 
from about 5 to 120. 

The Reynolds number is defined as an average value for the 
stacks for simplicity. It may be expressed as 

Re = p_ J 2H I u = 2^ 

tk 

where the dynamics viscosity of the fluid is /if and H f is the average 
fluid channel thickness. 

3. Results and discussion 

3.1. NTU as a function of Reynolds number 

The number of transfer units is given as a function of Reynolds 
number in Fig. 3 for the cases where k s = 10 W/(m*K), similar to 
gadolinium, and k s = 240 W/(m*I<), similar to aluminum, respec¬ 
tively. The trend is threefold. Firstly, the effective NTU is decreased 
as the standard deviation increases (due to the increased influence 
of the flow maldistribution). Secondly, at large Reynolds numbers 
the maldistribution has little overall impact on the NTU since the 
heat exchanger is overwhelmed by a flow rate too large to sustain 
sufficient heat transfer (see Eq. (6)). Thirdly, at very low Reynolds 
numbers the maldistribution is seen to have a decreased influence 
in the case with the larger thermal conductivity of the heat ex¬ 
changer solid (Fig. 3(b)), whereas for the case with a lower thermal 
conductivity (Fig. 3(a)) the impact of maldistribution is severe at a 
standard deviation of 5% or higher. 


These results show that heat transfer performance in non-uni¬ 
form channels is a complex function of geometry, materials and 
operating conditions. In cases where the heat exchanger is domi¬ 
nated by axial conduction rather than fluid convection, such as 
low Reynolds number or large solid thermal conductivity, maldistri¬ 
bution has little impact on performance and the heat exchanger can 
be modeled as a uniform stack of channels with acceptable accuracy. 
However, as the Reynolds number is increased for the same geome¬ 
try, the sensitivity to maldistribution becomes more significant. As 
an example, consider Fig. 3(b) and a standard deviation of 20%. When 
the Reynolds number is around 5 the influence of maldistribution is 
virtually negligible. However, if it is increased to a value of 20 the 
effective NTU is reduced to 25% of the nominal (ideal) NTU. This ef¬ 
fect may explain why heat exchangers in certain cases perform clo¬ 
sely to a uniform stack at a lower flow rate while exhibiting 
performance below expectations as the fluid flow increases. 

Considering lower thermal conductivity (Fig. 3(a)) it is seen that 
even at very low Reynolds numbers the flow maldistribution re¬ 
duces heat exchanger performance significantly. For the case shown 
here, even 5% standard deviation results in a significant deviation 
from the nominal. The major difference between Fig. 3(a) and (b) 
is thermal conduction between adjacent fluid channels, which is 
shown to improve overall heat transfer performance. 

3.2. NTU as a function of standard deviation 

The number of transfer units normalized to the ideal case is gi¬ 
ven as a function of the standard deviation in Fig. 4 for three differ¬ 
ent Reynolds numbers (5, 10 and 80, respectively) and two 
different thermal conductivities (10 and 240 W/(m*I<) K, Fig. 4(a) 
and (b), respectively). The results again show that the variation 
in plate spacing can greatly reduce heat transfer performance. It 
is, perhaps, of greater interest that the performance degradation 
is a strong function of both the thermal conductivity of the solid 
and Reynolds number. Performance of the heat exchanger with 
higher thermal conductivity is less dependent on the standard 
deviation and shows better performance for all cases modeled than 
that with lower thermal conductivity. Operating at the larger Rey¬ 
nolds number, however, is seen to result in a performance degra¬ 
dation that depends less weakly on thermal conductivity (as is 
also indicated in Fig. 3). This indicates that as the conduction be¬ 
tween channels is reduced compared to convection to the fluid (a 
lower value of NTU CT ) the heat transfer performance is reduced. 

At a large thermal conductivity and small Reynolds number (see 
Fig. 4(b)) the performance degradation from a high standard devi¬ 
ation is clearly negligible (also seen in Fig. 3(b)). In this case the 
thermal conduction is a dominating process in the heat exchanger, 
allowing temperature gradients caused by varying fluid velocities 
to be smoothed out. Considering the low thermal conductivity 
(Fig. 4(a)) case it is, however, clear that the performance is de¬ 
graded the most at the lowest Reynolds number. This is mostly 
due to the normalization of the NTU, as the lowest Reynolds num¬ 
bers will result in the highest ideal NTU values, which are most 
sensitive to variations in fluid flow. 

3.3. NTU as a function of thermal conductivity 

The NTU is given as a function of thermal conductivity at Rey¬ 
nolds numbers of 10 and 80, respectively, in Fig. 5. The plots are 
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(a) (b) 


Fig. 3. The effective number of transfer units, NTU st ack. (Eq. (17)) as a function of Reynolds number. The dashed line indicates ideal performance (homogeneous stack) and the 
legend otherwise refers to the standard deviations: (a) the solid thermal conductivity is 10 W/(m*K), (b) the solid thermal conductivity is 240 W/(m*K). 




Standard deviation [%] Standard deviation [%] 


(a) 


(b) 


Fig. 4. The effective number of transfer units (Eq. (17)) as a function of standard deviation normalized to the zero standard deviation case. The legend refers to the Reynolds 
number: (a) the thermal conductivity is 10 W/(m*I<), (b) the thermal conductivity is 240 W/(m*K). 


normalized to the ideal case for comparison. Again, the trend in 
heat transfer performance as a function of thermal conductivity 
does not fit a general form and is dependent on the specific heat 
exchanger. Considering the high Reynolds number case (Fig. 5(b)) 
it is seen that the thermal conductivity has relatively little influ¬ 
ence on the heat transfer performance and the trend is nearly lin¬ 
ear with thermal conductivity. At the higher flow rate, even large 
solid thermal conductivities do not improve heat exchanger perfor¬ 
mance markedly. However, the standard deviation is seen to have a 
significant impact. It may be concluded that in this case a simpler 
approach than the detailed, 2-dimensional stack model, could be 
applied for analyzing the heat exchanger performance. Each chan¬ 
nel may be treated individually and the final result would be a 
combination of single channel models. 

However, when considering the lower Reynolds number case of 
the same heat exchanger (Fig. 5(a)) a significant variation of the 
NTU as a function of k s is observed. At large and small values of 
the thermal conductivity the thermal cross talk has little influence 
(the graphs in Fig. 5(a) seem to be somewhat S-shaped). There 
does, however, exist a fairly broad range of intermediate values 
of the conductivity (from around 75 to, perhaps, 300 W/(m*K) in 


the case considered here) where the thermal cross talk is clearly 
significant and cannot be ignored or accurately treated with a sim¬ 
pler modeling approach. 

3.4. NTU as a function of cross talk 

The normalized NTU is plotted as a function of NTU CT for two 
plate spacing standard deviations in Fig. 6. For both standard devi¬ 
ation values, the normalized performance increases with increas¬ 
ing NTUcr- The curves tend to collapse onto themselves when 
plotted against NTU CT as opposed to solid thermal conductivity, 
suggesting that the NTU CT parameter is a sensible variable to con¬ 
sider when evaluating cross talk. Considering Fig. 6(a) when NTU CT 
is below a value of approximately 250, the heat transfer perfor¬ 
mance is a weak function of cross talk and the heat exchanger 
may be treated as a stack of independent channels that do not 
interact with one another, thus reducing model complexity. For 
intermediate values of NTU CT (250-1500) performance is highly 
dependent on interaction between flow channels and a fully de¬ 
tailed model of the stack is necessary to accurately predict heat ex¬ 
changer performance. For high values of NTU CT (above 1500), 
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k s [Wm'V 1 ] k s [Wm ’k" 1 ] 


(a) 


(b) 


Fig. 5. The effective number of transfer units (Eq. (17)) as a function of standard deviation normalized to the zero standard deviation case: (a) the Reynolds number is 10, (b) 
the Reynolds number is 80. The values are normalized to the zero standard deviation case for easier comparison. 




NTU 


CT 


(a) 


(b) 


Fig. 6. The effective NTU as a function of NTUcr (Eq. (20)) normalized to the zero standard deviation case at three different Reynolds numbers (indicated in the figure 
legends): (a) the standard deviation is 10%, (b) the standard deviation is 20%. 


variations in flow channel heights have a weak effect on perfor¬ 
mance and the heat exchanger can be treated as a uniform stack 
of channels. It should be noted that the values of NTUcr shown 
in Fig. 6(a) are not general and apply only to heat exchangers sim¬ 
ilar to the one modeled here. Fig. 6(b) shows a similar trend as that 
shown in Fig. 6(a) except the value of NTUcr necessaiy to be con¬ 
sidered a uniform stack increases to approximately 2500 because 
more cross talk between channels is necessary to smooth out tem¬ 
perature variations caused by the larger spacing variation. 

One important conclusion to draw from the study of Fig. 6 is 
that the level of modeling detail necessaiy to accurately character¬ 
ize a heat exchanger stack depends on a large number of system 
parameters, including the operating conditions. A heat exchanger 
that can practically be considered uniform when operating with 
a given solid material at a given flow rate may no longer agree with 
a uniform stack model if the flow rate is increased or it is con¬ 
structed with a solid with a lower thermal conductivity. Similarly, 
a given heat exchanger design may begin to deviate from a model 
that previously accurately predicted performance if the size is 


scaled down to a level where the manufacturing tolerance becomes 
a significantly larger fraction of the channel height. The term 
NTUcr can be used to give an idea of the significance of the thermal 
interaction between fluid channels and which type of modeling 
technique is most appropriate for a heat exchanger. 


4. Conclusion 

It was demonstrated that heat exchangers comprised of inho¬ 
mogeneous stacks of plates are subject to significant performance 
degradation over a large span of relevant parameters (Reynolds 
number, solid thermal conductivity and standard deviation of the 
parallel plate heat exchanger stacks) compared to uniform plate 
stacks. Only in cases of large thermal conductivity and low Rey¬ 
nolds number does the effect of flow maldistribution due to stack 
inhomogeneity become insignificant. 

It was shown that the heat exchanger performance may 
generally be divided into three dominating regimes. At low 














































































































































K.K. Nielsen et al. / International Journal of Heat and Mass Transfer 60 (2013) 432-439 


439 


thermal conductivity the stacks act as a bundle of individual chan¬ 
nels that are not communicating thermally (the cross talk is virtu¬ 
ally negligible). At elevated conductivities (somewhat depending 
on standard deviation) small changes in thermal conductivity 
may lead to significant changes in the performance reduction 
(caused by the flow maldistribution). Finally, at very large values 
of the thermal conductivity it was demonstrated that the effect 
of flow maldistribution is of minor importance and that heat 
exchangers operating in this regime may be treated as homoge¬ 
neous stacks. 

It is generally concluded that considering heat exchangers built 
of parallel plates (or similar geometries) may perform significantly 
less than predicted if using a model that assumes homogeneity in 
terms of perfect channels (and plates). Even at standard deviations 
of 5-10% the performance degradation due to flow maldistribution 
may be severe depending on the Reynolds number at which the 
heat exchanger is operated and the thermal conductivity of the so¬ 
lid. At larger solid thermal conductivities the effect of flow maldis¬ 
tribution may be reduced. 
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